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Abstract 

By treating the centers of solitons as point particles and studying their discrete dynamics, we demonstrate a new approach to the 
quantization of the soUton solutions of the sine-Gordon equation, one of the first model nonlinear field equations. In particular, 
we show that a linear superposition of the non-interacting shapes of two solitons offers a qualitative (and to a good approximation 
quantitative) description of the true two-soliton solution, provided that the trajectories of the centers of the superimposed solitons 
are considered unknown. Via variational calculus, we establish that the dynamics of the quasi-particles obey a pseudo-Newtonian 
law, which includes cross-mass terms. The successful identification of the governing equations of the (discrete) quasi-particles 
from the (continuous) field equation shows that the proposed approach provides a basis for the passage from the continuous to 
a discrete description of the field. 
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1. Introduction 

A sought-after property of model field equations is 
that they possess localized, permanent wave solutions 
that retain their identity upon interactions with one an- 
other. Solitons, which are solutions of fully -integrable 
equations, are an example of such waves. Unfortunately, 
integrability is not always a property of models that are 
of physical importance. Therefore, it is important to de- 
velop simple, albeit approximate, approaches to study- 
ing the dynamics of permanent waves in non-integrable 
systems, for which exact solutions are difficult to obtain. 
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The idea of identifying a localized solution of a non- 
linear wave equation with an elementary particle was 
first proposed by Perring & Skyrme They found a 
solution of (what is known today as) the sine-Gordon 
equation (SGE) 10 consisting of two interacting, lo- 
calized waves, which is an example of (what is nowa- 
days called) a two-soliton solution, several years before 
the notion of a soliton was introduced by Zabusky & 
Kruskal f3l. Through numerical simulation, the latter 
authors discovered that the localized traveling-wave so- 
lutions of the Korteweg-de Vries equation retain their 
shapes (identities) after they pass through each other (in- 
teract). Apparently, Zabusky & Kruskal were unaware 
of the work of Perring & Skyrme (|T| and arrived at the 
idea of identifying the nonlinear waves' dynamics with 
those of particles through experiment rather than con- 
jecture. 

Since their discovery, solitons have attracted an enor- 
mous amount of attention. Significant progress has been 
made in their mathematical description, and their appli- 
cations have been far-reaching |4 5 1. The so-called kink 
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solitary waves considered by Perring & Skyrme |[T] have 
been shown to be indeed soUtons |5 6|. Moreover, the 
relationship between the particle-like dynamics of the 
coherent structures that emerge in the solutions of non- 
linear wave equations and the field theories of particle 
physics are well-established in the literature li7<4.8i9J . 

From a mathematical point of view, we can elucidate 
the latter relationship by, somehow, reducing the "in- 
finitely complex" continuous description of the field to 
a "finitely complex" discrete description. In this Let- 
ter, we show how this can be achieved by studying the 
dynamics of quasi-particles . Our approach amounts to 
"degrading" the continuous description of the wave pro- 
file to a discrete description of the centers of coher- 
ent structures, assuming that the shapes of the coher- 
ent structures (for an integrable system, this would sim- 
ply be the solitons) are not significantly affected by the 
each other's presence. To this end, in this Letter, the 
term "quasi-particle" refers to these permanent, inde- 
structible and virtually non-deformable coherent struc- 
tures, whose centers can be treated as point particles 
with mass equal to some associated measure of iner- 
tia that we call the pseudomass. Furthermore, replacing 
a complicated continuous profile by a superposition of 
localized shapes can also be viewed as a coarse-grain 
description, in the sense that small deformations and 
wiggles are filtered out from the profile leaving merely 
the main structure (the "grains"). 

Hence, the coarse-grain description amounts to re- 
placing the solution of a nonlinear wave equation with 
a linear superposition of basis states (i.e., traveling- 
wave solutions in their undeformed, or non-interacting, 
state), then the trajectories of the centers of the super- 
imposed waves, which for a nonlinear equation do not 
follow the undisturbed (linear) trajectories, are treated 
as unknown. Restricting to the case of just two super- 
imposed waves, a discrete model for the trajectories is 
derived and solved numerically in this Letter This ap- 
proach gives a wave profile whose deviation from the 
analytical two-soliton solution (when available) is or- 
ders of magnitude smaller than the characteristic size of 
the solitons. Consequently, this paves the way to con- 
structing successive approximations that account for the 
higher-order, nonlinear interactions of solitons. 

Here, we note that our coarse-grain description is a 
special interpretation of the more general method of 
collective coordinates/variables, which has been put on 
solid theoretical ground [lOJ and become part of the 
textbooks on solitons Q. The latter approach made its 
debut in the study of resonances and collisions of soli- 
tary waves in the so-called 0* equation |11 12 1, a close 
relative of the SGE, and the study of two-soliton inter- 



actions in the SGE 0131141 . Moreover, the method of 
collective variables is just one type of variational ap- 
proximation, which is another approach to the analyti- 
cal study of nonlinear wave equations that has recently 
regained popularity [ 151161171 . Furthermore, it appears 
that Rice [181 was the first to realize that the collective- 
coordinate variational approximation provides a way 
of "distilling" the particle-like dynamics of nonlinear 
waves from the continuos (field) description, though 
Karpman & Solov'ev ifTsl had the foresight to use the 
term 'quasi-particle' in their discussion. 

Finally, we note that the nonlinear wave equation fea- 
tured herein — the sine-Gordon equation — continues 
to be of interest as a model field equation lfT9l[ . In ad- 
dition, various modifications of it have been considered 
in the literature. For example, in order to establish the 
effects of acceleration on the shape of the SGE's soli- 
tons, Fogel et al. |20| introduced a driving force into the 
SGE. However, this required also adding dissipation in 
the field equation in order to stabilize the evolution of 
the solitons, i.e., to ensure that they reach a steady ter- 
minal velocity 02OI21L Adding dissipation opens new 
horizons of investigation, and different physical mecha- 
nisms can be considered as progenitors of the dissipative 
force. It is well known that, in fluid mechanics, hnear 
dissipation of either viscous or Darcy type can balance 
the nonlinearity in the field equation, and allow stable 
localized waves to exist M22I23I . Moreover, in incom- 
pressible shallow-water flows, a viscous dissipation can 
allow for the existence of localized coherent structures 
with solitonic behavior ll24l[ . Nonetheless, one thing is 
certain: dissipation can alter the behavior of a nonlin- 
ear field equation dramatically. Therefore, in this Letter, 
we focus on the lossless SGE of Perring & Skyrme |[T] 
and show that the coarse-grain description of the field 
leads to the physical dynamics of locally -accelerating 
quasi-particles, without introducing a driving force or 
dissipation into the equation. In this respect, however, 
the SGE differs fundamentally from the modern loss- 
less nonlinear field theories of continuum mechanics 
(see, e.g., Ref. f25\ and those therein), since the "un- 
balanced" nonlinearity in the former allows for the cre- 
ation of localized coherent structures and does not lead 
to formation of singularities in finite time. 

2. The sine-Gordon equation and its soliton 
solutions 

For the purposes of this Letter, the SGE takes the 
following dimensionless (h = c = m = 1) form: 

u,t -Uxx = -sinu, (1) 



2 



where the subscripts denote partial differentiation. We 
have selected the SGE as our featuring example be- 
cause there are known analytical expressions for its two- 
soliton solutions. This allows us to show that the coarse- 
grain description is both an effective approximation tool 
and a new method for nonlinear-wave quantization. 

It is easy to show that the Lagrangian and the Hamil- 
tonian of the SGE read 



L,H= / juj^ (^^ul + I— cos u)(ix 



(2) 



respectively I12J5J, where the "— " sign in the right-hand 
side refers to L and the "+" sign to H. In addition, the 
wave momentum is defined |7| as 



P = 



u^Ut dx. 



(3) 



Then, the conservation of the energy and linear momen- 
tum require that AH /At = and AP/At = 0. 

Now, if we consider the moving frame ^ = x — vt, 
Eq. ([T]l reduces to an ODE, which has the following 
solution: 



u = (^((^) = 4arctan 



exp 



0< v< 1. 



(4) 

Notice that the latter is a one-soliton solution because 
2^[limt^ooM(x,f) — \\mx^-oou{x,t)] = 1 for all f < oo 
111. The kink, whose shape is termed a hydraulic jump 
in fluid mechanics, is a localized wave in the sense that 
its derivative 



: sech 



(5) 



is a function of finite energy (i.e., it is square-integrable). 

Another well-known analytical solution of the SGE is 
the two-soliton solution ||6l, which can take two distinct 
forms: 



= 4 arctan 



= 4 arctan 



cosh[(ei-e2)/2] ] 

fli2sinh[(0i + 02)/2]J' 

sinh[(ei-e2)/2] ] 
|fli2|cosh[(0i + 02)/2] J' 

where, for j G {1,2}, 

dj = yj{x^vjt-xj), Y^ = {\-x^)- 
a]^{\-Vj)l{\+v,), 
an = (fli -fl2)/(fli +fl2)- 



(6a) 



(6b) 



(7a) 
(7b) 
(7c) 



It is important to note that, due to the multivalued na- 
ture of the inverse tangent function, one has to be care- 
ful about the definition of m+. The inverse tangent in 
Eq. ( |6al ) should not be taken in the "principal value" 





Fig. 1. (Color online) Space-time plots of the soliton-.soliton (top) 
and the soliton-antisoliton (bottom) interaction. 

sense because, then, it leads to a discontinuous defini- 
tion of M+. But, if the quadrant of the argument of the 
inverse tangent is taken into account, then m+ is well- 
defined. 

Henceforth, the solutions m+ and are referred to 
as the soliton-soliton (S-S) and soliton-antisoliton (S- 
A) solutions, respectively. To justify this terminology, 
we present the spatial derivatives, which are the quanti- 
ties of importance in defining the corresponding quasi- 
particles, of the S-S and S-A solutions in Fig. [T] 

Now, when the two solitons are far from each other 
(as t +°o), it can be shown that Eq. (|6]) reduces to 

4 arctan [exp ( 0i + 51*= ) ] ± 4 arctan [exp ( 02 + ) ] , 

(8) 



where 5^ — —5^ — jlnjanl are the phase shifts Jl]. 



The above decomposition means that far from the point 
of collision of the two solitons, the profile is merely 
the linear superposition of two single solitons [whose 
shapes are identical to the one-sohton profile given by 
Eq. (|4|]. However, we cannot immediately express the 
two-soliton solution as a sum of two single solitons be- 
cause of the phase shifts, which are, in fact, the most 
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conspicuous manifestation of the nonlinearity as far as 
interacting solitons are concerned. Yet, the actual wave 
profile is approximated quite well by the mere sum of 
one-soliton profiles provided their centers follow trajec- 
tories that somehow account for the nonlinear interac- 
tion. The technical details and physical underpinning of 
this approximation are the goal of the this Letter. 

Here, we must note that Bowtell & Stuart f9] de- 
veloped a method, based on the poles of the Hamilto- 
nian density, that allowed them to compute the exact 
(ODEs for the) trajectories of the centers of the soli- 
tons. Their technique has not been applied to any other 
well-known nonlinear wave equations, which makes the 
SGE unique from the point of view of availability of 
reference solution for the trajectories of the centers to 
establish the practical applicability of the present ap- 
proach. However, we emphasize again that we do not 
intend to develop an exact method for computing the 
trajectories based on particular properties of the gov- 
erning nonlinear wave equation, rather we show how 
to approximate them based on more general principles. 
Thus, the method proposed below is applicable to inte- 
grable and non-integrable equations alike. 

3. From continuous to discrete dynamics 

3.1. One-soliton case 

To elucidate the concept of a quasi-particle we con- 
sider the propagation of a single solitary wave, whose 
shape is given by Eq. Q, centered at certain spatial po- 
sition x = X{t) (termed the particle trajectory), which 
may be unknown in advance. Under the latter assump- 
tions we have that 



M(x,f) = 0(^), £,=X-X{t), 



v{t)'^'x{t)- 



(9) 



Consequently, the time derivative of wave profile is 



-X{t) 



„2^3 



sech 



-Xit) 



(10) 

Now, the straightforward way to compute the discrete 
Lagrangian and Hamiltonian for the motion of the cen- 
ter of the soliton is to insert Eqs. ^ 8l ( fTO] ) into Eq. (|2| 
and to integrate over the infinite domain. However, in 
order not to obscure the main idea of the present work, 
we restrict ourselves to the non-relativistic (or classi- 
cal) approximation, i.e., X = v <^ 1 (small 
phase speed) and X = v = ^{v) (small phase acceler- 
ation). This means that the speed of the quasi-particle 
is much smaller than the characteristic speed of the 
equation (i.e., unity) and that the asymptotic order (of 



"smallness") of the acceleration is the same as that of 
the velocity. Under these assumptions the first term in 
Eq. ( [TO] i can be disregarded because it is of second order 
of "smallness" while the second term is of first order. 
In addition, any explicit dependence of the shape of the 
soliton on its phase speed is negligible in the classical 
limit. Finally, we note that the classical approximation 
is valid for speeds up to v ~ 0.2. 

Thus, under the above assumptions, we get the fol- 
lowing simple expressions for the profile and its partial 
derivatives: 

u = 4arctan[exp(^)], 1 — cosm — 2sech^^ , 
M;t = 2sech^, M, = — 2Xsech^. 

Then, by substituting the latter expressions into Eq. (|2]l, 
we get that the discrete Lagrangian L is 



(11) 



\MX^ 



, def 



sech2<^d^=8, (12) 



which is identical to the Lagrangian of a particle of mass 
M undergoing inertial motion. Similarly, we obtain, by 
substituting Eqs. ( [TT| ) into Eq. (|3]l, that the discrete wave 
momentum P is just 



-4Xsech2^d^ =MX, (13) 



which (again) is identical to the momentum of a (point) 
particle of mass M undergoing inertial motion. To dis- 
tinguish between the properties of the continuous and 
discrete descriptions, we call M the pseudomass and P 
the pseudomomentum. Hence, Eq. ( [T3| ) shows that the 
wave momentum of a single solitary wave is the mo- 
mentum of its "center," if the latter is considered as a 
point particle endowed with mass M = 8. Of course, 
this should be of no surprise if the solutions to the SGE 
are to model elementary particles, which obey the prin- 
ciple of wave-particle duality. 

So far we have used the term "center of a soliton" 
without defining it. Unfortunately, there is no general 
agreement on the matter, and the center of a coherent 
structure (e.g., a soliton) has to be defined on a case-by- 
case basis. A definition is given by Bergman et al. f26l 
but, as Kaup |27| shows, depending upon the definition 
being used, one may obtain that the quasi-particle's dy- 
namics are Newtonian |20| or non-Newtonian f2Tl — 
clearly, this is quite a discrepancy. In this respect, the 
closest work to ours is that of Rice ITSl . who derives 
the discrete Hamiltonian and Lagrangian for a single 
quasi-particle whose support can be considered an in- 
ternal parameter of the problem. He comes to the con- 
clusion that the dynamics of a deforming quasi-particle 
are, generally, non-Newtonian, as do we. For the pur- 
poses of this Letter, it suffices to define the center of 
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a coherent structure u as the point where Ux attains its 
maximum (recall Fig. [TJ, and that is precisely where 
the corresponding quasi-particle should be located. Note 
that, for the SGE solitons, this also corresponds to the 
inflection point of u{-,t), which is always located at x* 
such that Uxx{x*,t) — 0. 



3.2. Two-soliton case 

Continuing on to the two-soliton case, we first note 
that one can formally represent any two-soliton (two- 
wave, if the equation is non-integrable) solution of the 
field equations as the combination of the undeformed 
(non-interacting) shapes of two solitary waves <I>i, <I>2 
and an arbitrary disturbance term «l>i2 due to the inter- 
action of the latter two: 



u{x,t)^^l[x~Xi{t)]+^2[x-X2{t)] 

+ ^l2[x^Xx{t),X-X2{t) 



(14) 



In fact, one can think of <I>i2 as the difference between 
the full two-soliton solution given by Eq. (|6| and the lin- 
ear superposition of two single soliton given by Eq. (|8]l, 
with properly selected trajectories that are unknown in 
advance. 

The two first terms of Eq. ( [I4j i give a linear superpo- 
sition that is bound to be a poor quantitative description 
of interaction of two nonlinear waves when they overlap 
significantly. Indeed, it is well-known that, e.g., for the 
Korteweg-de Vries equation, the maximal amplitude of 
a two-soliton solution is much smaller than the sum of 
the maximal amplitudes of the individual solitons |5|. 
But, a linear superposition along two undisturbed tra- 
jectories Xi{t) = Xi{Q) + v,r would give a maximal am- 
plitude equal to the sum of the two maximal amplitudes. 

The main idea of the present work is that if the 
trajectories X,(f) are selected in a specific manner, 
then the predominant part of the nonlinear effects, 
represented by <I>i2 in Eq. ( [T4] l, in the two-soliton so- 
lution can be recovered by the linear superposition 
alone, i.e., Eq. ([T4| with <I>i2 neglected. That is to say, 
for properly chosen trajectories, one has max|<I>i2| < 
emax{max|<I>i|,max|<I>2|} for some small e(> 0). It 
is not know in advance how small e is, but one can 
estimate it a posteriori and vindicate the approxima- 
tion. We do so by comparison to the exact solution of 
the problem, but one can also judge the quality of the 
approximation, without any knowledge of the exact so- 
lution, by the novel method of Kaup & Vogel \7E\ . The 
latter will be necessary when the governing equation is 
non-integrable and the former approach fails. 



Now, we show how the trajectories of the quasi- 
particles [as quantified by the functions X\ (?) and X2(f )] 
can be calculated. First, recall that we are working in 
the classical limit so we can neglect the dependence of 
a soliton's shape on its phase speed and use the unde- 
formed (non-interacting) shape of the soliton in the cal- 

def 

culation. Thus, after setting <I>i = <I>2 = <t>, neglecting 
the term <I>i2 in Eq. ( [T4] i and acknowledging Eq. (|8j, we 
arrive at the following simple ensemble of two waves 
as the coarse-grain approximation: 

u^{x,t)^^[x~Xy{t)]±^[x^X2{t)l 
4>(Q =4arctan[exp(C)], 

where the "+" and " sign follow the convention set 
forth in Eq. (|6|. Then, the first partial derivatives of the 
wave profile are just 

uf{x,t) = -X,^'[x-X^{t)] TX2^'[x~X2{t)], (I6a) 

u^{x,t)= 0'[x-Xi(f)]± 4>'[x-X2(f)]. (16b) 

Henceforth, we leave the time dependence of X\ and 
X2 implicit to simplify the notation. Consequently, the 
kinetic energy term of the Lagrangian, for both m+ and 
(denoted simply by m), reads 



(15) 



/: 



Awf djc= \MiiXl + \M22Xl 



+ Mi2(X2-Xi)XiX2, (17) 

where the pseudomasses M/y (/, y S { 1 , 2}) are given by 



12(X2-Xi) = ±/ 4>'(x-Xi)<I>'(x-X2)dx 



= ±4 y sech(x — Xi)sech(x — X2)dx 
= ±8 (X2 - Xi ) csch(X2 - ) , 

(18a) 

\^\^i)fAx = 4 |^°°sech2^,d^,- - 8, 

(18b) 

and ^; = x — X\. 

Similarly, we can compute the interaction potential 
(i.e., the remaining terms of the Lagrangian). First, we 
note that 

jM^ + 1 — cosM = 4{sech^(x — Xi) + sech^(x — X2)± 

[cosh(2;c - Xi - X2 ) - 1 ] sech^ (jc - Xi ) sech^ (x - X2 ) } . 

(19) 

Then, by integrating the latter [Eq. ( [T9| l] over spac^^ 
we find that the potential energy of interaction is 



^ These integrals are rather difficult, thus we refer the reader to 
the appendices of Ferguson & Willis 1141 and Sugiyama 111 111 for 
explicit formulae and hints. 
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f/ (z) = 1 6 ± 4 csch^ ( iz) sech ( ^z) (sinhz - z) , (20) 

def 

and we have set z = X2—X1 for convenience. Conse- 
quently, combining the above computations, we obtain 
the discrete Lagrangian by substituting the coarse-grain 
approximation, given by Eq. ( [T5] l, into Eq. (|2|, namely 

lL=\MnXf + Mn{z)XiX2 + \M22Xi~U{z) (21) 

Thus, we can use the Euler-Lagrange equations for 
the extremization of the discrete Lagrangian to find the 
governing system of equations for the centers of the 
solitons (i.e., the quasi-particles), namely 

^ (MnXi) + ^ [Mniz)X2] +M'i2(z)^i^2 - U'{z), 

(22a) 

^ [Mn{z)Xi] + ^ (M22X2) -M'i2(z)^i^2 = -U'{z). 

(22b) 

Finally, recalling that M„- are constants, X? <^ 1 for / S 
{1,2}, and using the chain rule for differentiation, we 
can recast the above system into 

MnXi+Mniz)X2= U'{z), (23a) 
Mn{z)Xi+M22X2 = -U'{z). (23b) 

The most conspicuous trait of the above equation of 
motion is that it is not Newtonian per se. It reflects 
the fact that the inertia of each particle is affected by 
the presence of another (accelerating) particle in its 
vicinit}p] Therefore, it is more appropriate to refer to 
such pseudo-Newtonian dynamics as Machean. Here, 
by "Machean dynamics" we mean those that satisfy a 
local version of Mach 's principle, which is the hypoth- 
esis that the inertia of a body depends not only upon its 
mass but also upon the quantity, distribution and rela- 
tive motion of matter in its viscinity |30|. In Eqs. ( |23] l 
this effect manifests itself through the terms whose co- 
efficient is Mi2(z). In fact, the crossmass Mi2(z) gives 
rise to an interaction between the quasi-particles simi- 
lar to the interaction due to the potential U{z), so the 
discrete kinetic and potential energies (i.e., those of the 
quasi-particles) are not in one-to-one correspondence 
with their continuous counterparts (i.e., those of the 
solitons). 

Now, it is convenient to render Eqs. ( |23] l into "New- 
tonian form" by resolving them with respect to the sec- 
ond time derivatives of each quasi-particle, namely 



' This should not alarm the reader, as it is now generally accepted 
that Newton's concept of inertia and his second law of motion are 
merely the starting point for modern theories (see, e.g., Ref. 1291 
and those therein), not the end-all. 



Xi= f/'(z)E(z)(M22 + Mi2) = G(z), (24a) 
X2 = -f/'(z)S(z)(Mn +M12) = -G(z), (24b) 
where 

= (MiiM22-M?2)"' (64-64z2csch2z)"' 

(25) 

for the case when Mn = M22. For this same situation, 
the "mathematical" (or effective) force per unit mass 
takes the form 

G{z) = ±[z — 2sinhz— (sinhz — 2z) cosh z] 

X 2csch4(iz) sech2(iz) /(8 TSzcschz), (26) 

where the "+" and "— " sign still refer to the soliton- 
soliton and soliton-antisoliton case, respectively. Here, 
it should be noted that the asymptotic behavior of the 
right-hand side, as z ^ 0, is G(z) ^ z^' for the S-S case 
and G(z) ~ z for the S-A case. The singularity (i.e., 
the infinite repulsion) in the S-S case is precisely what 
prevents the trajectories of the quasi-particles (solitons) 
from crossing each other 

Finally, we note that Ferguson & Willis lfT4l previ- 
ously performed a number of calculations, in the context 
of Josephson junctions, similar to the ones presented 
above. Nonetheless, our results are novel because we 
propose a very different understanding of the results ob- 
tained by the collective-variable method, one that has 
fundamental physical meaning [embodied by Eqs. (|23]l]. 

4. Numerical results 

4. 1 . Quasi-particle dynamics 

In this section, we provide numerical results, based 
on the coarse-grain description of the soliton-soliton 
and soliton-antisoliton interactions, and benchmark 
them against the exact results in f9l. Though a quan- 
titative comparison cannot be made because no scales 
are provided on the plots in ||9 |, we can still perform a 
qualitative comparison. 

It is clear that the system given by Eqs. (|24]l can be 
reduced to a single differential equation governing the 
evolution of the distance between the quasi-particles, 
i.e., X2—X\. Then, upon obtaining the solution for 
X2—X1, one can go back and solve Eqs. ( [24] l for each 
individual trajectory. For the sake of simplicity, and 
convenience in comparing the results to the analytical 
two-soliton solution, we consider here the case of two 
quasi-particles whose initial phase speeds and positions 
are equal in magnitude but differ in sign. Consequently, 
we have 
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Xi(f) = y(f), X2(f) = -y(r), (27a) 

Xi(0) = y(0) = -xo, X2{0) = -Y{0)= xo, (27b) 
Xi(0) = y(0)= vo, X2(0) = -7(0) = -vo, (27c) 

and, clearly, X2{t) —Xi (f ) = ~2Y{t). Now, we can com- 
bine Eqs. ( p4j i and ( p7] i to obtain the following initial- 
value the problem for Y{t): 

Y = G{-2Y), Y{0)^-xo, F(0) = vo. (28) 

The analytic expression for the effective force (per unit 
mass) G{z) is known and given by Eq. ( |26| l, so we 
can use any standard numerical ODE integration algo- 
rithm, e.g., Mathematica's NDSolve, to solve the ini- 
tial value problem for Y{t) given in Eq. ( |28| ). Naturally, 
when the function Y{t) is obtained numerically, we can 
recover the trajectories of the two quasi-particles. 

For the initial phase speed, we select a value that 
falls within the classical limit, namely vq = 0.1. Respec- 
tively, we choose an initial position of xq — 6, which 
verifies the assumption that (initially) the interaction 
between the two solitons is insignificant, i.e., they are 
weU-separated. 

Fig. |2] presents the results for the S-S (top) and S- 
A (bottom) solutions. Recall that we defined the cen- 
ters of the solitons, i.e., the points where has a local 
maximum (see Fig. [T|i, as the locations of their corre- 
sponding quasi-particles. So, the difference between the 
straight lines in Fig. |2j which represent the trajectories 
followed by two non-interacting quasi-particles (soli- 
tons), i.e., Xj — ^XQ ± VQt, and the actual trajectories 
at the time when the quasi-particles (solitons) return to 
their initial positions gives their phase shifts. 

For the S-S case, according to the numerical solution, 
the time needed for the quasi-particles to return to their 
starting position [i.e., t* ^0 such that F(f*) = F(0)] is 
t* 74.0355. On the other hand, the time needed for 
the latter quasi-particles to reach the same spatial po- 
sitions would be f* = 120 in the absence of interaction 
between the two quasi-particles. Hence, the interact- 
ing quasi-particles return to their initial positions At = 
120 — 74.0355 = 45.9645 dimensionless time units ear- 
lier than their non-interacting counterparts. Thus, if the 
non-interacting quasi-particles were to return to their 
initial positions at the same time as the interacting ones, 
the former would have to have been separated by a dis- 
tance voAf = 0. 1 X 45.9645 = 4.59645 smaller than the 
actual distance of 12 dimensionless length units (see 
Fig.[2]l. Therefore, half of this number, 2.29823, is pre- 
cisely the phase shift of each quasi-particle. At the same 
time, substituting the values of the initial phase speed 
(i.e., Vo = 0.1) in Eqs. Q & ([8|, we can compute the 
exact phase shift of each of the solitons in the full two- 
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Fig. 2. (Color online) The thick solid lines represent the trajectories 
of the interacting quasi-particles (as calculated from coarse-grain 
description) are shown for the S-S case (top) and the S-A case (bot- 
tom). The dashed lines represent the trajectories of two non-inter- 
acting quasi-particles (solitons). And. the thin solid lines represent 
the location of the inflection points of the two-soliton solution. 

soliton solution (for both the S-S and S-A case): 

5i = -52 = 2.30259. (29) 

Clearly, the phase shift we obtained from our numerical 
computation is in very good quantitative agreement with 
the exact one — the difference is merely 0.2%. 

For the S-A case, the solution for the trajectory is 
shown in Fig. |2] (bottom). The time needed for the 
quasi-particles to swap their positions after their inter- 
action is t* = 75.9069, which is Af = 120 - 75.9069 = 
44.0931 units earlier than their non-interacting counter- 
parts. Therefore, the phase shift of each quasi-particle 
is 2.20466, which gives a relative error between the 
approximate and exact phase shifts of 4%. Indeed, in 
this case, the approximation is supposed to be gener- 
ally worse than in the S-S case, since in the S-A case 
the solitons overlap significantly and their velocities get 
large. This wiU become clear in the next section. 

Thus, the coarse-grain description gives a good, 
quantitative account of the trajectories of the centers of 
solitons (i.e., the location of the corresponding quasi- 
particles). This provides further support for the claim 
that solitons are indeed particle-like objects and, to a 
good approximation, they can be considered as point 
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particles obeying pseudo-Newtonian dynamics. 

4.2. Soliton dynamics 

Apart from the significance of the present work 
as tool for studying the discrete dynamics of quasi- 
particles corresponding to solitons, our results also 
show that the coarse-grain description is a practical way 
of constructing approximate wave-train solutions (two 
waves in this Letter but, eventually, A^) of nonlinear 
wave equations. This is of great importance for non- 
integrable systems where analytical expressions for 
wave-train solutions are not available. Here, we outline 
the quantitative limits of this kind of approximation 
by comparing our results to the analytical two-soliton 
solution of the SGE. 

To this end, we compute the difference between the 
analytical solution, given in Eqs. (|6]l, and the coarse- 
grain approximation, given by Eq. ([TSj, presented in the 
previous section. Fig. [3] (top) shows the results for the 
S-S case, where the maximal pointwise error is found 
to be w 0.028, which corresponds to a relative error of 
0.028/(4 ;r) w 0.23%. In addition, we note that, for the 
S-S case, the maximum of the phase speed occurs at the 
initial instant of time, hence the speeds of the solitons 
(and quasi-particles) remain classical for all times under 
consideration. 

The S-A case is depicted in Fig. |3] (bottom) and the 
maximal pointwise error for this case is « 0.73, which 
corresponds to a relative error of 0.73 / {2n) « 12%. The 
reason the error for the S-A case is an order of magni- 
tude larger than the errors for the S-S case is that the 
assumption of classical speeds is no longer justified. 
In particular, the attraction between the quasi-particles 
increases their velocities (i.e., the phase speeds of the 
corresponding solitons) significantly. As a result, we 
have max, v(f) =Xi (37.9534) = 0.675316, which is no 
longer small in comparison to unity. Nonetheless, even 
with the relativistic effects ignored, the error in the 
approximation (away from the point where the quasi- 
particles annihilate each other) is small enough so that 
the coarse-grain description is in good agreement with 
the exact two-soliton solution. 

Here, we note that, for the purposes of these compar- 
isons, we have eliminated any phase error by evaluat- 
ing the two-soliton solution at f — tJ*, where tJ* is the 
phase error committed by the approximation. Thus, the 
instant of time when the quasi-particles are closest (S-S 
case) or pass through (S-A case) each other predicted 
by the coarse-grain description is the same as the one 
given by the analytical solution. We found that iJ+ = 




Fig. 3. (Color online) Space-time plots of the absolute value of 
the difference between the exact two-soliton solution and the linear 
superposition of two single solitons following the calculated trajec- 
tories for the S-S case (top) and S-A case (bottom). 

-0.0660218 for the S-S case and = 0.858073 for the 
S-A case. Finally, to make sure the quasi-particles and 
the solitons start at the same place, we had to take, in the 
analytical solution given in Eqs. (|6]), xi2 = =f3 -70838 
for the S-S case and xi^ = T3.70954 for the S-A case. 

5. Conclusions 

An approximate approach to the construction of two- 
soliton solutions of nonlinear wave equations is dis- 
cussed. It consist of considering the linear superposition 
of the non-interacting (or inertially-propagating) shapes 
of two single solitons following trajectories that are un- 
known in advance. On the basis of the variational formu- 
lation of the field problem, a discrete Lagrangian is de- 
rived that governs the trajectories of the corresponding 
quasi-particles (i.e., the point-particles endowed with 
certain mass located at the position of the center of a 
soliton). 

The approach was validated using the celebrated sine- 
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Gordon equation, which has analytical two-soliton so- 
lutions and the dynamics of the corresponding quasi- 
particles have been investigate. By numerically integrat- 
ing the ODEs that govern the trajectories of the quasi- 
particles, we constructed an approximate two-soliton 
profile, showing that the quasi-particles obey a pseudo- 
Newtonian, more precisely Machean, law. Then, the 
quasi-particle dynamics and the coarse-grain solution 
were compared to the exact dynamics [9J and the ana- 
lytical two-soliton solution, respectively. In both cases, 
excellent quantitative and qualitative agreement was ob- 
tained. Furthermore, the present results can easily be 
extended to the relativistic case by including another 
collective-coordinate that governs the width of the wave 

mm. 

Finally, we emphasize that the collective-coordi- 
nate/variable (or variational approximation) type of 
approach discussed herein is an effective technique for 
studying the two-soliton solutions of equations that are 
not fully-integrable. Therefore, our interpretation of the 
method — i.e., the idea of a coarse-grain description — 
can shed light on the passage form the continuous to a 
discrete description of the dynamics (i.e., the quantiza- 
tion) of nonlinear dispersive waves, their relation to the 
elementary particles and the concept of wave-particle 
duality for more than just the "classical" (integrable) 
equations of soliton theory. 
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